Determination of the 5d6s^Di state lifetime and blackbody radiation clock shift in Yb 
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The Stark shift of the ytterbium optical clock transition due to room temperature blackbody 
radiation is dominated by a static Stark effect, which was recently measured to high accuracy 
[J. A. Sherman et ai, Phys. Rev. Lett. 108, 153002 (2012)]. However, room temperature operation 
of the clock at 10~^* inaccuracy requires a dynamic correction to this static approximation. This 
dynamic correction largely depends on a single electric dipole matrix element for which theoreti- 
cally and experimentally derived values disagree significantly. We determine this important matrix 
element by two independent methods, which yield consistent values. Along with precise radiative 
lifetimes of 6s6p^Pi and 5d6s^Di, we report the clock's blackbody radiation shift to 0.05% precision. 

PACS numbers: 06.30.Ft,32.60.+i,32.70.Cs 



Alkaline-earth-like atoms, such as Yb [T], Sr [2H1], and 
Hg [S] feature intrinsically narrow ^Sq -H- ^Pq optical 
transitions capable of serving as stable and accurate fre- 
quency references [6] when cooled and held in an optical 
lattice trapping potential [3 [S] . Accurate knowledge of 
clock transition frequencies advances timekeeping tech- 
nology and enables new tests of physics [SHTT] . 

Atomic frequency references are defined by an ideal 
system: atoms at rest in a null-field, zero-temperature 
environment |12j . If a physical realization deviates from 
these ideals, researchers must account for corrections to 
the measured transition frequency and, importantly, un- 
certainty present in these corrections. Here, we explore 
the dominant ytterbium clock correction [H I13| due to 
room-temperature blackbody radiation (BBR). 

The polarizing effect of BBR largely mimics that of a 
static electric field due to the low frequency nature of 
BBR relative to optical transitions involving clock states 
(see Fig. [T]). Writing the BBR clock frequency shift [14] 



Aj/bbr = - 



1 Aa{0) 



{E^)t[1 + Vcioc^iT)], (1) 



highlights its similarity to a static Stark shift, where 

Aa(0) = ae(0) - ag(0) = 145.726(3) a.u., (2) 

is the differential static polarizability between clock 
states \g) = |6s^ ^Sq) and |e) = |6s6p^Po), now known to 
high accuracy [T3] (a.u. = atomic units [IS]). {E^)t ~ 
(8.3193 V/cm)2(T/300 K)"^ is the time-averaged electric 
field intensity of BBR at absolute temperature T [IB] . 
A small dynamic correction 77ciock(300 K) < 0.02 ac- 
counts for the frequency dependence of Aa(aj). Over 
90% of ?7ciock depends on coupling between 6s6p ^Pq and 
neighboring 5d6s ^Di (sec [14 and supplemental material 
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FIG. 1. (color online) Most of the room-temperature BBR 
energy spectrum (thick red line) is far infrared of transitions 
involving the Yb clock states. The states' polarizabilities 
[see Eq. ([3|], largely constant over much of the BBR spec- 
trum, are balanced at the 'magic' trapping frequency uj* . 



(SM)). But, critically, a measurement [TT and recent pre- 
cise calculation 18J of the electric dipole matrix element 

V= \{6s6p^Po\\T>\\5d6s^Di)l 

differ by ^ 3 standard deviations. 

In this work, we independently determine V in or- 
der to accurately compute ?/ciock and the clock correc- 
tion Ai^BBR- We present two distinct approaches re- 
sulting in good agreement. First, we describe a semi- 
empirical technique which combines existing polarizabil- 
ity data with atomic theory to constrain V. Then we de- 
scribe a measurement of the 5d6s '^Di radiative lifetime 
in trapped Yb. Since 5d6s ^Di decays predominantly to 
the 6s6p ^Pj manifold, V is readily extracted. Finally, we 
discuss current accuracy limitations imposed by Az/bbr- 

Method I: Semi- empirical technique — Accurately mea- 
sured experimental parameters, such as the differential 
static polarizability [Eq. ([2|], also depend on electric 
dipole coupling between 6s6p '^Pq and 5d6s '^Di , and sub- 
sequently can be used to constrain the value of T) |19| . 
The electric dipole polarizability of clock state n due to 
radiation at an off- resonant frequency uj/2n reads 



a„(w) 



1 2 



E 



\{n'\mn) 



(3) 



2 



where (n'||D||n) is a reduced electric dipole matrix ele- 
ment and a;„'„/27r = (Wn' — Wn) /h is the corresponding 
transition frequency. As w —> 0, we recover the familiar 
static polarizability expression. 

The so-called 'magic' trapping frequency uj* , which 
balances the polarizabilities of the clock states 6J, 



Aa{uj*) = ae(w*) — ag{u!*) ~ 0, 



(4) 



has been measured to high accuracy in Yb [U [20] . Equa- 
tions ^ and Q may be combined to yield 



Aa(0) + bAa{cu*) = 145.726(3) a.u.. 



(5) 



where b is arbitrary but may be chosen to to our benefit. 
In the linear combination 

a„(0)+6«„(c.*) = ^2 ^ |(„'||D||n)|' — 



;i3 , 

X ( 1 + 6 "^"''^ 



up' , — UJ*^ 
n'n 



(6) 



the term in parenthesis serves as a 'scale factor' rela- 
tive to each transition's static polarizability contribu- 
tion. For instance, for a choice 6 = — 1, this scale fac- 
tor tends to zero for uj^'n » UJ* because these transi- 
tions contribute nearly identically to both polarizabil- 
ities a„(0) and a„(w*) [see Eq. ([3)]. We find advan- 
tage in choosing a value 6 — 1 such that contributions 
from certain low-lying transitions are suppressed in the 
linear combination Aa(0) -I- bAa{uj*), along with con- 
tributions from the higher-lying transitions. In Table 



IV 



we present contributions to Eq. Q from the lowest- 
lying transitions in each clock state for both 6 = and 
b — —0.75. In each case, we write the contribution 
from 6s6p^Pq 5d6s^D^ in terms of unspecified ma- 
trix element T>. Other contributions are derived from 
experimental lifetimes in Refs. pTI - E5] . Transitions to 
the closely spaced states (4/^'^)5d6s6p (|, |) ^ and ■^P^ 



('mixed states' in Table IV) are exceptions; these con- 
tributions were estimated with a CI-I-MBPT calculation 
similar to Ref. [TH]. It is evident from Table IV how 



the choice of b affects the relative importance of cer- 
tain transitions. For example, while the 'mixed states' 
contribute sizably to the differential static polarizability 
Aa(0), their contribution to Aa(0) — 0.75Aa(a;*) is neg- 
ligible. Moreover, contributions from higher-lying transi- 
tions not explicitly shown in the table — which contribute 
at the ~ 10% level for both state polarizabilities [21] — 
are also largely suppressed with the choice b = —0.75. 
Specifically, the 'scale factor' in Eq. ([6| is nearly zero for 
the lowest of these transitions (for which w„'„ ~ 2 a;*), 
rising to just 0.25 for the highest-lying, least important 
transitions. 

Tallying contributions from all transitions, we find 



Aa(0) - 0.75Aa(tJ* 



26.8X>2 -64(8) +0(6), (7) 



in atomic units. Here the first two terms on the r.h.s. ac- 



TABLE I. Contributions to the static polarizability q,i (0) and 
linear combination a„{0) — 0.75Q:n(a;*) [refer to Eq. m] for 
the lowest-lying transitions from the clock states (a.u. 



Qn(0) 



Q„(0) - 0.75a„(aj* 



6s6p ^Pi 
6s6p ^Pi 



5d6s ''Di 
6s 7s ^Si 
6s6d^Di 
6s8s '■'Si 
mixed states 



clock state n — 6s Sq 

2 -1 

100 -4 

21 1 

clock state n — 6s6p ^Pq 

20.W^ 26.8©^ 

37 -65 

22 -3 
2 
39 



the uncertainty is dominated by that of the matrix el- 
ement {6s7s^Si\\D\\6s6p^Po) f23 . The additional term 
0(6) accounts for contributions from all higher-lying tran- 

We ascribe an 



sitions not given explicitly in Table IV 



uncertainty to this term based on experimental upper 
limits to the polarizabilities of the two clock states [24], 
along with theoretical input from Ref. [TH] and present 
CI-HMBPT calculations. Equating the r.h.s. of Q to 
experimental result ([5| gives V = 2.80(7) a.u. We com- 
pare this result with other determinations and new data 
below. 

Method II: lifetime measurement — Alternatively, mea- 
surement of the 5d6s ^Di radiative lifetime yields T> 
since = {3TTeohc^Co){2J' + l)/(wgT„), where J' = 1; 
ujq/2tt « 2.1587 x lO" Hz and Co = 0.64(1) are the radi- 
ated frequency and branching ratio to ^Pqj respectively. 
Co is accurately computable because LS'-coupling remains 
vahd 121]. In the cascade 5d6s^Di 6s6p^Pi 65^ i^'o 
[see Fig. [2]^a)], atoms emit a 556 nm photon during the 
second decay which is technically easier to detect than 
the first radiated (infrared) photon [17]. Other states 
populated by the decay {^Pq and '^P2) are long-lived. If 
atoms are instantaneously excited at time to to ^Di, flu- 
orescence from ^Pi follows a double exponential [ITl [26] , 



count for contributions from all transitions in Table IV 



y{t) =Axe{t~ to) [e-f*-*")/^" - e-'*-*")/^" 1 + yo, (8) 

where is the radiative lifetime of "^Pi (r^ > r^), A is 
a scaling factor and i/q accounts for stray detected light. 
The Heaviside unit-step Q{t — to) models rapid atom ex- 
citation at <o- Decay branching ratios affect only the 
normalization of Eq. Q , not its time dependence \TT^ . 

We describe the cooling and confinement of iVat ~ 10'* 
atoms of ^^^Yb in a one-dimensional optical lattice else- 
where pj. As depicted in Fig.[2jb), a resonant '7r-pulse' of 
578 nm light [2H] coherently transfers atoms from ^Sq to 
the long-lived '^Pq state. Then, a brief (te ~ 25 ns) reso- 
nant pulse of 1388 nm light excites more than half of these 
atoms to '^Di. An event counter accumulates the arrival 
times of radiated 556 nm photons into 5 ns bins. Heat- 
ing due to photon scattering, background gas collisions. 
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FIG. 2. (color online) We measure the 5d6s ^Di radiative life- 
time via cascade decay, a) Key lifetimes are labeled Ta and 
Ti,. Double-arrowed lines indicate laser transitions. Wavy 
lines are important decay channels, b) A representative pulse 
timing diagram, c) Typical fluorescence data, fit (red line) to 
Eq. ([sjl. An inset highlights the signal and fit near an excita- 
tion time to — 523.3(6) ns. d) Observed lifetimes vary with 
atomic density, p. Lines are linear regression fits. Logarithmic 
scaling emphasizes data at low p with negligible interaction 
effects. Error bars represent standard uncertainties obtained 
from nonlinear fits, and uncertainty in p estimates. 



and accumulation in '^P2 limit the number of excitations 
per loading cycle to about 200. Though we estimate a 
modest light collection/detection efficiency (^ 0.1 %), we 
typically observe iVat x (1.5 x 10~^) green photons per 
excitation. Biases due to photon 'pile-up' in counter bins 
are negligible. Between 10^-10^ excitations are sufficient 
to obtain satisfactory decay profiles [e.g., Fig.[2]jc)]. 

We fit fluorescence signals to Eq. (8]) with a 
statistically-weighted Levenberg-Marquardt routine. 
Though covariance in A, Ta, and T{, can be significant, 
simulations establish that fitting biases become negli- 
gible with sufficient count totals. In large data sets, 
event counter technical noise synchronous with a timing 
oscillator overwhelms the signal shot noise. Re-binning 
data into 20 ns chunks removes much of this noise, but 
does not significantly alter the results. A maximum 
likelihood method yielded statistically similar fits. 

Atomic interactions, such as collective emission ('su- 
perradiance', 'subradiance') or radiation trapping may 
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FIG. 3. (color online) The present measurement (*) and semi- 
empirical result (f ) (open red circles) are compared to other 
determinations. Measurements a-j are from, respectively, 
[29H36] ■ |17) . [37] • We note that the small error bar on h ac- 
counts only for statistical uncertainty. Green points k-m are 
calculations [38] , [18] , [39] . When necessary, we infer lifetimes 
from reported matrix elements or natural linewidths. We as- 
sign an error bar to the Ta prediction of / using the authors' 
estimate of uncertainty in a highly correlated polarizability. 



TABLE II. The uncertainty in lifetimes Ta and rj, due to 
atomic interactions is largely statistical since we report ex- 
trapolations to zero atomic density. Covariant fit biases and 
distortions due to Zeeman oscillations are estimated from 
Monte-Carlo simulations. Uncertainties due to 1388 nm pulse 
duration (r^) and stray light are statistically limited. 





u{Ta) (ns) 


u{Tb) (ns) 


Atomic interactions 


4.3 


3.3 


Fit biases 


0.9 


1.5 


Zeeman quantum beats 


3.0 


3.0 


1388 nm finite excitation 


3.3 


4.3 


1388 nm stray light 


3.4 


4.6 


Event counter timing 


0.2 


0.4 


Total (quadrature sum) 


7.1 


7.4 



influence radiative decay. We probed these eflFects by 
varying the atomic density p undergoing decay [4^ . Re- 
sults [Fig. [2jd)] indicate non-negligible shortening of Ta 
at high p. In the limit of slow dipole dephasing [IT], col- 
lective emission shortens an observed decay lifetime as 
T = To(l + pLA^/4)~-^, where tq is the single atom value, 
L is the length of a pencil-shaped atomic cloud, and A 
is the radiated wavelength [35]. For p — 10^ cm~'^ and 
i = 0.1 mm, the modification in this simple model is 
about 5%. We see the opposite effect in ^Pi, an increase 
of Tb at high p. We explored both effects by varying 
578 nm and 1388 nm excitation pulse areas, altering the 
relative populations of -^Sq, ^Pq, and ^Di, but observed 
no substantial change in the systematic effect; further 
investigation is warranted. Nevertheless, the effects of in- 
teractions on observed lifetimes are made negligible over 
the lowest decade of examined densities [see Fig. [2jd)]. 
Quenching due to cold-collisions and lattice scattering is 
negligible. 

Potential systematic effects arise from finite 1388 nm 
pulse duration Te, and spurious excitation due to poor ex- 
tinction. We varied Tg from 25 ns to 90 ns and observed 
no significant systematic effect on fitted state lifetimes. 



4 



TABLE in. The present results, expressed as reduced matrix 
elements (a.u.), are compared with selected literature values. 





V 


|{6s6p3Pi||D||6s2i5'o)| 


Experiment 


2.77(4) 


0.542(2) 


Semi-empirical 


2.80(7) 




Experiment [TtJ 


2.58(10) 


0.547(16) 


Calculation [38] 


2.61(10) 


0.54(8) 


Calculation [TS] 


2.91(7) 


0.587 


Calculation [39) 


2.58(23) 


0.41(1) 



We configured the apparatus for fast actuation and high 
(60 dB) extinction of the 1388 nm light. A single-mode 
fiber-coupled mechanical shutter (500 /iS rise-time) pre- 
cedes a fiber-coupled acousto-optic modulator (AOM). 
A second (free-space) AOM, driven with a tuned delay, 
aids in pulse shaping and extinction. A beam dump, iris, 
and baffles minimize the influence of scattered light. The 
final AOM imposes a 160 MHz frequency shift, detuning 
scattered light from resonance and increasing effective 
extinction. With up to 1 mW of deflected 1388 nm light 
focused to a 30 ^m spot (an intensity I ~ 7x 10^ W/m^), 
we attain a Rabi frequency $7 cx I'D exceeding 1 GHz. 
We observed no significant changes to Tq and Tb when 
varying an optical attenuator over 4 dB to test for de- 
pendence on Q and stray light intensity. 

We systematically excited to both hyperflne compo- 
nents ^Di{F' = ^,^), which are split by 3.07(7) GHz. We 
observe no hyperfine quantum beats .43j due to the large 
splitting and selective laser excitation. We observed no 
Zeeman oscillations |17j or significantly different results 
when an applied magnetic field B varied from 0.01 mT 
to 0.1 mT. 1388 nm excitation light propagated along 
the lattice axis and was polarized perpendicular to B. 
We detected 556 nm photons ~ 45° from the lattice axis 
with largely polarization insensitive optics. We observed 
a slight dependence on the number of excitations per 
loading cycle but believe this effect is due to atomic inter- 
actions since increased scattering reduces the atom num- 
ber. Varying the lattice laser intensity over 50% yielded 
no significant change in observed lifetimes. 

We report the radiative lifetimes Tq = (329.3 ± 7.1) ns 
and Tft = (866.1 ± 7.4) ns. Table [TT| enumerates measure- 
ment uncertainties. We compare our results to existing 
measurements and calculations in Fig. [3] Our result for 
Tb agrees with many prior measurements. Our result for 
Ta lies between the only other measurement [T7| and a 



recent calculation [TH]. Our semi-empirical method ex- 
hibits good agreement with the measurement. Table HI 
lists the results as inferred matrix elements. 

BBR clock shift—The BBR Stark shift to the clock 
frequency is found from the expression 



AfBBR 



1 



2he, 



Jo 



(9) 



where uy(cj) is the BBR spectral energy density cor- 
responding to temperature T, given by Planck's law. 
A static approximation neglecting the slight frequency 
dependence of Aa(a;) over the BBR spectrum (refer 
to Fig. [l]) is formally obtained by making the sub- 
stitution /S.a{u}) — >■ Aq:(0) in Eq. An improved 
approximation takes into account the lowest-order fre- 
quency dependence of the polarizability arising from the 
low-lying 6s6p'^Pg — ^dGs^D^ transition: Aa(w) — >■ 
Aa(0) + {2/Zh){'D'^ /ujDuj'^. Integrating over w analyt- 
ically, we interpret the additional shift as f7ciock(r) ~ 



0.017 (. 



from Eq. (^1|, where 



63 (hua)^ Aq(0) ~ V300 K,* 

ks is Boltzmann's constant. A more thorouglT account 
of small contributions from all other transitions, includ- 
ing the ^5*0 state and next-order terms (oc T^) yields 



r/ciock(T) = 0.0173(5) 



T 



300 K 



0.0006 



/ T 



V300 K 



We omit higher order terms (oc T^,T*,...) which are 
negligible at T < 300 K. We provide more details of this 
evaluation, including multipolar effects [Tl], in the SM. 

Conclusion — Assuming an ideal BBR environment at 
300 K, we use the present results to calculate Az^bbr = 
— 1.2774(6) Hz. The present determination of jydock sets 
an uncertainty limit for Ai^bbr at 1.1 x 10~^^. In prac- 
tice, uncertainty in Az^bbr also arises from imprecise 
knowledge of the thermal environment. In an existing 
apparatus, we estimate an effective temperature uncer- 
tainty of 1 K due to non-uniformity, corresponding to 
a fractional clock uncertainty of 3.3 x 10~^^ [T]. The 
present results therefore motivate further efforts to con- 
trol the thermal environment of the clock chamber |44j . 
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I. TABULATED RESULTS FOR ATOMIC FACTORS 



Table IV presents transition frequencies and dipole matrix elements for transitions from the clock states, as well 
as corresponding contributions to static polarizability a„(0), polarizability at the magic frequency a„(w*), and the 
linear combination a„(0) — 0.75a„(cij*) for both clock states. Also provided are contributions to the atomic properties 

(2) 

an and Pn, which are discussed in the following section. 

In the upper portion of Table [TV] explicit values are given where experimental data is available, with uncertainties 
being derived from the corresponding references. For remaining transitions, contributions to the static polarizability 
are given as unknown parameters Xi,2,3- Corresponding contributions to the other atomic properties are a written as 
a range (a, 6)xi.2,3, where the range (a, 6) is determined purely by the atomic spectrum and magic frequency. This 
helps illustrate the relative importance of the unknown contributions in each case. For example, while X2 denotes 
the combined contribution to static polarizablity ae(0) due to the 'mixed states', the associated contribution to 
ae(0) — 0.75ae(cj*) is comparatively suppressed, necessarily being between —0.033x2 and 0.004x2- 



In the middle portion of Table IV contributions to the differential properties [e.g., Aa{0) = ae{0) — ag(0)] are 
tallied. Here A('main states') incorporates all values given explicitly in the upper portion of the table. A('mixed 
states') is estimated by CI+MBPT calculations similar to those described in Ref. (TS], whereas A('all others') is 
estimated with CI+MBPT calculations together with additional theoretical input from Refs. [181 [2l]. Theoretical 
uncertainty for A('mixed states') and A('all others') is difficult to assess; the numbers given in the table represent 
reasonable estimates of this uncertainty. We reiterate that methods developed in the main text minimize the influence 
of these contributions. 



In Table IV contributions from the 6s6p "^Pq 5d6s ^D^ transition are derived from our most accurate determination 
(see main text) of the matrix element V = \ {6s6p^P„\\'D\\5d6s^Di) \ = 2.77(4). In Table |vj we compare these results 
with contributions derived from other determinations of this matrix element. 



TABLE IV. Transition frequencies, matrix elements, and contributions to various atomic properties of interest for the Yb clock 
states. The magic frequency is uj* — 0.06000 [T][20|. {a,b)x denotes a range (ox, &x)- All values are in atomic units. 



n' 




|(n'llDlln)| 




an(0) 




a„(0)-0.75a„(a;*) 




Pn/10^ 












clock state 


n = 6s^ % 






6s6p ^Pi 


0.08198 


0.542 ± O.QOi 


a 


2.39 ±0.02 


5.14 ±0.04 


-1.47± 0.01 


0.355 ±0.003 


-0.219 ± 0.002 


6s6p ^Pi 


0.11422 


4.148 ±0.002 




100.40 ±0.09 


138.7 ±0.1 


-3.601 ±0.003 


7.696 ± 0.007 


-0.2760 ± 0.0003 


(4/»)5d6s^ (1,1) 


^ 0.13148 


2.03 ±0.0-^ 


r 


20.8 ±0.8 


26 ± 1 


1.10 ± 0.04 


1.20 ±0.05 


0.063 ±0.002 


all others 




Xi 


(1.00,1.14)xi 


(0.15,0.25)xi 


(0.00,0.03)xi 


(0.000, 0.005)xi 












clock state 


n = 6s6p^Po 






5d6s ^Di 


0.03281 


2.77± 0.0_f 




156 ±4 


-67 ±2 


206 ±5 


145 ±4 


192 ±5 




0.07020 


2.0±0.]R 




37 ±4 


136 ± 16 


-65 ±8 


7.4 ±0.9 


-13 ±2 




0.10261 


1.82±0.0# 




21.6 ±0.7 


33 ±1 


-3.02 ±0.09 


2.05 ±0.06 


-0.287 ±0.009 


6s8s^Si 


0.11084 


0.64 ± 0.0-^ ' 


F 


2.4 ±0.3 


3.4 ±0.4 


-0.15 ±0.02 


0.20 ±0.02 


-0.012 ± 0.001 


mixed states 








X2 


(1.33,1.38)X2 


(-0.033, 0.004)x2 


(0.07,0.08)X2 


(-0.0025, 0.0003)x2 


all others 








X3 


(1.00,1.31)X3 


(0.02,0.25)X3 


(0.00,0.07)X3 


(0.000, 0.005)x3 


A (main) 








93 ±6 


-64± 16 


142 ± 10 


146 ±4 


179 ±5 


A (mixed states) 








39± 13 


52± 17 


0.0 ±0.5 


2.7 ±0.9 


0.1 ±0.1 


A (all others) 








16± 16 


20 ±20 


0±6 


1 ± 1 


0.08 ±0.08 


A (total) 








148 ± 21 


8 ±31 


142 ± 11 


149 ±4 


179 ±5 


A(expt.) 








145.726 ± 0.003 





145.726 ± 0.003 







expt., present 
^ Ref. 12] 

weighted mean from four values compiled in Ref. 
Ref. [23] 
" Ref. \m 
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TABLE V. Comparison of the 6s6p "^Pq — > 5d6s "^Di contribution to the atomic properties in Table IV with different deter- 
minations of the dipole matrix element. Here e specifies the excited clock state 6s6p^Pq. All values are in atomic units. 





(6s6p^PollD|l5d6s^L>i) 


ae{0) 


ae(a;*) 


ae(0) -0.75ae(c^*) 


ai^VlO' 


Pe/10^ 


unspecified 


V 


20.320^ 


-8.6660^ 


26.82D^ 


18.882?^ 


24.92D^ 


expt., present 


2.77 ±0.04 


156 ±4 


-67 ±2 


206 ±5 


145 ±4 


192 ±5 


semi-emp., present 


2.80 ±0.07 


159 ±8 


-68 ±3 


210± 11 


148 ±7 


195 ± 10 


expt., Ref [n] 


2.58 ±0.10 


135 ± 10 


-58 ±4 


179 ± 14 


126 ± 10 


166 ± 13 


theor., Ref [H] 


2.91 ±0.07 


172 ±8 


-73 ±4 


227 ± 11 


160 ±8 


211 ± 10 


theor., Ref [39] 


2.58 ±0.23 


135 ± 24 


-58± 10 


179 ± 32 


126 ± 22 


166 ± 30 



II. BBR CLOCK SHIFT 

Atomic units are employed throughout this and the following section. The usual definitions, e = = h = 47reo = 1, 
are supplemented with the additional definition fc^ = 1, where fc^ is Boltzmann's constant. We also define a reference 
temperature Tp equivalent to 300 Kelvin and having a value Tq = 9.50 x 10^'* in our system of units. We base our 
unit system on SI electromagnetic expressions, with the Bohr magneton being given by /is — eTi/2me — 1/2. The 
speed of light c is used in favor of the fine structure constant a in expressions to follow to avoid notational confusion; 
c — ~ 137 in atomic units. 

The energy shift to clock state n due to electric dipole coupling with thermal radiation reads 



Ai?„ = — 27r / u{u]) an{uj) duj , 
Jo 



(10) 



where u{uj) is the spectral energy density at (angular) frequency uj and ctn(^) is the frequency-dependent polarizability, 

(jJji'n 



a„{uj) 



I |(n'||D||n)|' 



(11) 



The Cauchy principal value is implicitly taken for the integral of Eq. (10), as well as for the integral of Eq. (|13[) to 



follow [47j . For BBR, the spectral energy density u{uj) may be written in terms of temperature T using Planck's law, 



u{u!) 



1 



7r2c3 e"/"?^ - 1 ' 



and it follows that the energy shift may be recast as 



where F(y) is the function introduced by Farley and Wing |47| . 

2 f°° 



1 1 

h 

y+x y-x 



X 

- 1' 



(12) 



(13) 



The function F{y) is displayed in Figure |4j The atomic transition frequencies all satisfy ujn'n S> Tq, and it follows 
that for room temperature (or below), the factor F{ujn'n/T) appearing in Eq. (12) may be well-approximated with 
the leading terms of the asymptotic expansion F(y) = 47r^/452/ ± 327r^/189y^ ± 327r''/45?;^ + . . . [T|1|37]- The energy 
shift may then be decomposed into respective terms, 



2^4 
15c3 



v(0) 



63c3 



,(2) 



where we have introduced the frequency-independent atomic factors aii , 



15c3 

(fc) 



,(4) 



|(n'||D||n) 



k+1 



(14) 



(15) 
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FIG. 4. (color online) Function F{y) (solid) and it's asymptotic form truncated at the leading three terms (dashed). For Yb, 
all transition frequencies satisfy uj^'n/To > 34. We find the three-term asymptotic expansion of F{y) to be accurate to better 
than 4 parts in f 0^ for y > 34. 



Noting the equivalence between ai? and the static polarizabiUty, al^"^ = a„(0), we identify the leading term in 
Eq. ( 14) with the 'static approximation' of the energy shift. To connect with Eq. (1) of main text, we write A^'bbr = 
{AEe — Ai?g) /2tt with the static contribution factored out: 



(2) 



407r^T(f Aa 
21 Aa(0) 



4 4Aa(4) fT 
"^"^°A^l7i; 



+ 



(16) 



where Aaf*^' 



a 



(k) 



is the differential atomic factor taken between the excited (e) and ground (g) clock states. 



The correction factor 77ciock(r) can then be equated to the terms in square brackets succeeding the leading 1. 



III. EVALUATING 77ciock(T) 



With Aa^"^ known accurately from recent e xpe rimental measurement [13^, the leading term of ryciock is determined 
by the differential atomic factor Aa*^^^. Table IV tabulates contributions to Aa^^\ yielding a final result 

Aa(2) = 1.49(4) X 10^ 

We see that the 6s6p^Pq — 5d6s^Di transition dominates. While other contributions are smaller in comparison, they 
are nevertheless non-negligible at our accuracy. 

As mentioned previously, theoretical uncertainty is difficult to assess for contributions from 'mixed states' and 'all 
others'. For this reason, we also consider an alternative extraction of Aa*^^' by defining a property 

p„ ^ + [a„(0) - a„(^*)] = ^ ^ |(n'||D||n)|' ^ " 0.75 V^^.) 

From this definition, we get an expression for Aa*-^^: 

Aa(2) = Ap_ _ [Aq,(o) - Aa(a;*)] Ap - 30.3659(6) x 10^ 

where we have utilized accurately known experimental parameters Aq!(0), /S.a{ui*)^ and uj* in the last equali ty. Since 

gives 



IV 



the last term is known precisely, uncertainty in Aa^^^ is then commensurate with uncertainty in Ap. Table 
contributions to Ap, and we find 'mixed states' and 'all others' to give negligible contribution in this case. With Ap 
as given in Table |IV[ we then obtain 

Aa(2) ^ 1.48(5) X 10^ 

which is in agreement with our result given above. 

Using Aa^^) — 1.49(4) x 10^ we find the leading term to the dynamic correction factor to be: 

^^^- 0173(5) 
21 Aa(o) " "■"1^^^^-' 
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Contributions to q;„ scale as l/u'^^. This explains why 6s6p^Po ~^ ^dQa^D^ — whose transition frequency is less 
than half of all other transition frequencies — has an increased relative importance in Aa^^^ compared to Aa^'^^. For 
Aa'-^-', we find that Qs&p'^Pq — > bd&s '^D^ completely dominates, and we obtain for the next highest term of the dynamic 
correction factor: 

S^'r^^^ = 0.00059(1). 
Higher-order terms are found to be negligible. Finally, putting these results together, we find 

?7ciock(r) = 0.0173(5) (^Ij^^ +^-^^^^{^Y^ . 
where uncertainty in the second term may be neglected. 



A. Magnetic dipole and higher-multipolar couplings 

Thus far we have confined our attention to dominant electric dipole [El) coupling to the BBR field. Additionally, 
the atom couples to the BBR field via magnetic dipole (Ml) and higher multipolar (iJ2, M2, . . . ) interactions. Porsev 
and Derevianko [2] argued that for room temperature BBR, the Ml coupling could potentially lead to fractional 
frequency shifts on the level of 10~^* in optical lattice clocks. Higher multipolar couplings were shown to be suppressed 
below this level, and we neglect them here. 



he 



The Ml BBR shift is analogous to the El shift; it is given by Eq. (10) with the substitution and Q!„(aj) — > /3„(w)/c^ 
with j3n(uj) being the frequency-dependent magnetic polarizability. f3n(!-o) is defined analogously to Q!„(w) [Eq. (11 
but with the magnetic dipole operator /x replacing the electric dipole operator D. The additional factor c? in t 
substitution accounts for the different magnitudes of electric and magnetic fields in the BBR spectrum. 

Due to parity selection rules (/x is an even-parity operator, whereas D is an odd-parity operator) the Af 1-allowed 
transitions differ from the _El-allowed transitions. In particular, the ^Pq clock state has a low-frequency Ml transition 
to the neighboring ^P^ state of the S8-nie fine structure niRuifold. For this trcinsition, ujn'71 

/To sa 3.4, and it follows 

that the asymptotic expansion of F{y) is not appropriate (see Fig|4|. In this case, Eq. ( [l2| (with D — > /x/c) should 
be used directly. 

The Ml shift may be estimated by assuming the non-relativistic limit and absence of configuration mixing between 
states. In the non-relativistic limit, the magnetic dipole operator is given by the expression fi = — /iB(L -I- 2S), where 
L and S are the total orbital and spin angular momenta of the electrons, respectively. In the absence of configuration 
mixing, it follows that the only non-vanishing Ml matrix element involving either clock state is the one connecting 
the ^Pq clock state to the nearby ^Pi state. This matrix element may then be evaluated analytically, with the result 
|(nsnp'^PQ||/x||nsnp'^P]^)| — \/2/^s- Within these approximations, the room temperature Ml BBR shift to the ^Pq 
clock level is found to be 

where Wfs is the fine structure interval between the '^Pq and "^Pi states, and the factor F (wfs/To) ~ 0.28 for Yb. 
Interpreted as an additional contribution to f7ciock(2o), we find: 

.n^(T^o)-lxlO-. 



We therefore conclude that Ml coupling to the BBR field is negligible. 



